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Abstract. In this article we derive differential recursion relations for the Laguerre functions on the 
cone Q of positive definite real matrices. The highest weight representations of the group Sp(n,R) play 
a fundamental role. Each such representation acts on a Hilbert space of holomorphic functions on the 
tube domain Q + iSym(n, R). We then use the Laplace transform to carry the Lie algebra action over to 
L^(n, dfj.,^). The differential recursion relations result by restricting to a distinguished three dimensional 
subalgebra, which is isomorphic to s[(2,R). 



Introduction 

The theory of special functions has its origins in the late eighteenth and early nineteenth centuries 
when it was seen that the algebraic, exponential, and trigonometric functions (and their inverses) were 
not adequate to express results to differential equations that arose in the context of some important 
physical problems. New functions arose to which we have associated names like Bessel, Hermite, Jacobi, 
Laguerre, and Legendre. Then there are the Gamma, Beta, Hypergeometric and many other families of 
special functions. By the latter half of the 19*^ century these same functions arose in different contexts 
and their name 'special' began to take on greater meaning. Their functional properties were explored 
and included functional relations, differential and difference recursion relations, orthogonality relations, 
integral relations, and others. 

The preface of Vilenkin's book ^21 notes that the connection between special functions and group 
representations was first discovered by E. Cartan in the early part of the 20*^ century. By the time 
Vilenkin's book appeared in the 1960's that interplay had been well established. The texts by Miller, 
Vilenkin, and Vilenkin and Klimyk |10[ 1171 [T5j . for example, well document the general philosophy. In 
short, group representation theory made it possible to express the classical special functions as matrix 
entries of a representation and to unify many of the disparate relationships mentioned above. The 
representation can then be used to derive differential equations and differential recursion relations for 
those functions. 

In 1964 Simon Gindikin published his paper: 'Analysis on Homogeneous Domains', cf. jBj. This 
important paper developed special functions as part of analysis on homogeneous convex cones and built 
upon the earlier work of C. L. Siegel ^HI on the cone of positive definite matrices. The Siegel integral of 
the first and second kind generalize to become the Beta and Gamma functions for the cone, respectively. 
Generalized hypergeometric functions are extended to homogeneous cones. Many important differential 
properties also extend. 

Around the same time M. Koecher [7| |S] began to develop his analysis on symmetric cones and the 
complex tube domains associated with them. Jordan algebras proved to be a decisive tool for framing 
and obtaining many important fundamental results. The outstanding text by Faraut and Koranyi [S] 
documents this interaction (see also its extensive bibliography). Nevertheless, the representation theory 
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of Hermitian groups, which are naturally associated with tube domains, is not used in any outstanding 
way. 

In a series of papers ^ El El El the second and third authors (with Genkai Zhang in the first 
two referenced articles) use the representation theory of Hermitian groups in a decisive way to obtain 
differential and difference recursion relations on series of special functions. In the context of bounded 
symmetric domains the relevant special functions are generalized Meixner polynomials and in the context 
of tube domains over a symmetric cone the relevant special functions are Laguerre functions. These special 
functions exist in distinguished L^-spaces, which are unitarily isomorphic to Hilbert spaces of holomorphic 
functions on either a bounded symmetric domain or a tube domain. The well-known representation theory 
that exists there then transfers to the corresponding L^-space to produce differential and difference 
relations that exist among the special functions. One cannot downplay the essential role that Jordan 
algebras play in establishing and expressing many of the fundamental results obtained about orthogonal 
families of special functions defined on symmetric cones. Nevertheless, the theory of highest weight 
representations add fundamental new results not otherwise easily obtained. 

In this present paper we will continue the themes outlined in the above mentioned papers for the 
Laguerre functions defined of the cone of positive definite real symmetric matrices. The underlying group 
is Sp(n,R) and its representation theory establishes new differential recursion relations that Laguerre 
functions satisfy. The case n = 1 reduces to the classical Laguerre functions defined on R+. Briefly, the 
classical Laguerre polynomials are defined by the formula 

o.)^ 1^1^ (:)(-.)' 

and the Laguerre functions are defined by 

The classical differential recursion relations are then expressed by the following three formulas: 

(1) xD^ + vD- x)t„,{x) = -{2m + j^)C(2;) 

(2) xD"^ + {2x + v)D + {t + i^)K,ix) = -2m{v + m - 

(3) xD^ - {2x - u)D + {t- = -2C+i(a;)- 

It is these three formulas that we generalize via the representation theory of Sp(n, R) to Laguerre functions 
defined on the cone of positive definite real symmetric matrices. (Some definitions of LJ^ include a factor 
of 1/n!. This is the case in pP but not |21. The inclusion of this factor changes the differential recursion 
relations slightly.) 

This article is organized as follows: In the first section we introduce some standard Jordan algebra 
notation. In particular, we introduce the Laguerre functions and polynomials. Even though this material 
and most of the material in Sections 2 and 3 hold in general for simple Euclidean Jordan algebras, we 
specialize to the case J = Sym(n, M), the Jordan algebra of symmetric real n x n- matrices. In Section 
2 we introduce the tube domain T(S1) = fl + iSym(n,R), where fi is the open self dual cone of positive 
definite matrices. We also discuss the structure of the group Sp(n,R) and its Lie algebra sp(n, R). Some 
important subalgebras of sp(n, C) are introduced. This structure is later used to construct the differential 
operators that give rise to the differential equations satisfied by the Laguerre functions. 

Section 3 is devoted to the discussion of the highest weight representations {iT^,H„{T{n))). We also 
introduce the Laplace transform as a special case of the restriction principle introduced in [13 . In Section 
4 we describe how the Lie algebra, 5p(ri, C) acts on 'H,^{T{n)). In particular. Proposition 14.21 gives an 
explicit formula for the action of the derived representation for each of the three subalgebras 6c, p^, and 
p~, whose direct sum is sp(n, C). It should be noted, however, that not all of this information is needed 
to establish the differential recursion relations for the Laguerre functions. In fact, only the action of three 
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elements are needed. The action of sp(n,C) on L^{Q,d^,y) is described in Sectional The result is the 
following theorem: 

Theorem 15. 21 For / G Lf,{^) a, smooth vector we have: 

(1) X^{X)f{x) = tr[{bx + {ax - xa - vh)^ - xVV^]f{x), X ^ !) ^ 

(2) K{X)J{x) = tv[{va + ax+ {ax + xa + va)\I + xVaV]f{x), X = (^"^^ e p+ 

(3) A^(X)/(a;) = tr[(i/a - ax + {ax + xa - ua)W - xWaW]f{x), X = G p- 

Here we use V to denote the gradient, is the complexification of the Lie subalgebra u(n) C sp(n, M), 
and p* are certain Abelian subalgebras of sp(n,C) on which ic acts. We explain the main ideas for the 
special case of Sp(l,]R) ~ SL(2,R) and present the lengthy proof of the theorem in the Appendix. 

Specializing the above results to the elements 

where 1 stands for the n x n identity matrix, using properties of highest weight representations, and 
employing Lemma 5.5 of 2 we get our main result: 

Theorem 16. 31 The Laguerre functions are related by the following differential recursion relations: 

(1) tr(-a;VV - i/V + x)t^{x) = {nv + 2\xn\)t^{x). 

(2) tr(xVV + {vl + 2.x) V + {vl + xW^{x) = -2E;=i (m'^ej) ^"^^^ - 1 + ^ - (j " l))C-ej(-^) 

(3) tr(-a:VV + {-vl + 2x)V + {vl - x))C(a;) = 2 E;=i c™(j)C+e, {x). 



1. The Jordan Algebra of Real Symmetric Matrices 

In this section we introduce the Jordan algebra J of real symmetric matrices. We then discuss the space 
of L-invariant polynomial functions on J and the F-function associated to the cone of symmetric positive 
matrices. Finally we introduce the generalized Laguerre functions and polynomials. 

1.1. The Jordan Algebra J = Sym(n,IR). We denote by J the vector space of all real symmetric 
n X n matrices. The multiplication x o y = ■^{xy + yx) and the inner product {x\y) = ti{xy) turn J 
into a real Euclidean simple Jordan algebra. The determinant and trace functions for J are the usual 
determinant and trace of an n x n matrix and will be denoted det and tr, respectively. Observe that 
dim J -.^ d — !li!Ltll_ Lg^ fi denote the interior of the cone of squares: {x^ | x G j}. Then Q is the set of 
all positive definite matrices in J. Let 

H{n) = {g e GL( J) \gn = n} 

and let H be the connected component of the identity of 7?($7). Then H can be identified with GL(n, M)+ 
(where + indicates positive determinant) acting on 17 by the formula 

g ■ X — gxg^, g E H, x E il. 

This action is transitive and, since Q is self-dual, it follows that f2 is a symmetric cone. Let L be the 
stability subgroup of the identity e Eft. Then L = SO(n, R) and 



(1.1) 



n ~ H/L 
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Let Ei^i be the diagonal n x n matrix with 1 in the (i, i)-position and zeros elsewhere. Then 
{Ei^i, . . . , En.n) is a Jordan frame for J. Let J*-*^-* be the +l-eigenspace of the idempotent Ei^i + - ■ ■ + Ek^k 
acting on J by multiplication. Each j''''^ is a Jordan subalgebra and we have 

c c • • • c J(") - J. 

If detfc is the determinant function for Jk and Pk is orthogonal projection of J onto J*^*^^ then the function 
A/c(x) — detfc Pk{x) is the usual fc**^ principal minor for an n x n symmetric matrix; it is homogeneous of 
degree k. In particular A(a;) := A„(a;) = det(a;). Note also that 

(1.2) A(/i • x) = det(/i)2A(a;) , Vh e H . 

For m — (toi, . . . , m„) e C" we write m > if each uii is a nonnegative integer and nii > m2 > 
■ • ■ > "T-n ^ 0. We let A = {m | m > 0}. For each m G A define 

A A mi — m2 A ™2— "13 a "'n- 1 — '"71 a m„ 

— ^1 ^2 ■■■^ri-1 ■ 

These are the generalized power functions. It is not hard to see that the degree of Am is |m| := 
TOi + • • • + m„. Observe that each generalized power function extends to a holomorphic polynomial on 
Jc = Sym(n, C) in a unique way. 



1.2. The L-invariant Polynomials. For each m S A we define an L-invariant polynomial, V'm, on Jc 

by 



ipmiz) = J A^{lz)dl, z e Jc 
where dl is normalized Haar measure on L. A well known theorem of Schmid (c.f. |15j ) gives 

Lemma 1.1. IfV{Jc) is the space of all polynomial functions on Jc and V{Jc)^ denotes the space of 
L-invariant polynomials then {V'm | m > 0} is a basis of V{Jc)^ ■ Furthermore, if VkiJc)^ denotes the 
space of L-invariant polynomials of degree less than or equal k then {V'm | |m| < k} is a basis ofPkiJc)^- 

This lemma implies among other things that ipm{e + x) is a linear combination of V'n, |n| < |m| . 

The generalized binomial coefficients, (^^^ j ^^'^ defined by the equation: 

V'm(e + x)= (^)V'n(a;). 

n|<|m| ^ ^ 

1.3. The Generalized Gamma Function. For m e C" we define Am(a;), x ^ ft, hy the same formula 
given above for m G A. The generalized Gamma function is defined by 



rf,(m) = / e-*'-"Am(x)A(x)^^ dx. 
Jn 



Conditions for convergence of this integral are given in the proposition below. If A is a real number we 
will associate the multi-index (A, • • • , A) and denote it by A as well. The context of use should not cause 
confusion. Thus we define 

_ ro(A + m) 
^^^^ - To (A) • 

For later reference we note the following facts about the generalized Gamma function: 
Proposition 1.2. Let the notation be as above. Then the following holds: 
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(1) If m ~ (mi,TO2, . . . ,77in) G C" then the integral defining the generalized Gamma function con- 
verges if Re(mj) > (j — 1)^, for j — 1, . . . ,n, and in this case 

r^,(m) = (2^)^ n - - 1) 2), 

where T is the usual Gamma function. In particular it follows, that the T -function has a mero- 
morphic continuation to all of C"' . 

(2) // efc is an n-vector with 1 in the k*^ position and 's elsewhere then the following holds for all 
m e C"; 

(a) r7Sr^ = ™^-i-(^-4- 

lj2(m-efc) 2 
rj2(m + efe) 1 

Proof. Part (2) follows immediately from (1) and part (1) is Theorem 7.1.1 of [1]. □ 

1.4. The Generalized Laguerre Functions and Polynomials. Let v > and m € A. The general- 
ized Laguerre polynomial is defined (cf. 5 p. 242) by the formula 

^m(2^) = (t^)m XI 73-'/'n(-^)' ^ ^ J 

|n|<|m| ^ / ^" 

and the generalized Laguerre function is defined by 

CW = e-*--L^(2x). 

Remark 1.3. In the case n = 1, i.e. in the case G ~ 5^(1,1$.) = SL(2,M), the generalized Laguerre 
polynomials and functions defined above are precisely the classical Laguerre polynomials and functions 
defined on M+. We refer to P for the discussion of that case. 

The determinant Det(/i) oi h E H acting on J is 

Det(/i) = det(/i)"+i . 

It follows from H1.2|l that the measure 

d^oix) = A(x) ~ dx 

is iJ-invariant. Here dx is the Lebesgue measure on J. More generally, we set d^^{x) = A.{xY ^dx 
and define 

Ll{n) = L^{n,d^i„). 

We observe that by p. 2(1 it follows that H acts unitarily on Ll{Vl) be the formula 

K{h)f{x)^<:\ei{hYf{h^-x). 

Theorem 1.4 The set {P^ | m > 0} forms a complete orthogonal system in i^(il)^, the Hilhert 

space of L-invariant functions in Lj(ri). 

In it was shown, that the classical diff'erential recursion relations and diff'erential equations for 
the Laguerre functions on K+ follows from the representation theory of Sp(l,M) = SL(2,IR). In |3] this 
was generalized to the space of complex Hermitian matrices. It is a goal of this article to extend this 
result to the generalized Laguerre functions defined of the cone of symmetric matrices. This indicates 
nicely what the more general results should be. Here we use heavily the structure of Sp(n, M) and its Lie 
algebra, but the proof of the general results should be more in the line of Jordan algebras. 
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2. The Tube Domain T{n), the Group Sp(n,K) and its Lie Algebra 

In this section we introduce the tube domain T{Cl) — Cl + iSyni(n, M) and discuss the action of the group 
Sp(n, M) on this domain. We then discuss some important Lie subalgebras of sp(n, C), the complexification 
of the Lie algebra of g = sp(n,R). These subalgebras will show up again in Section 0] where we compute 
their action on Hilbert spaces of holomorphic functions on T{n) introduced in the next section. We then 
use that information to construct the Laguerre differential operators. 

2.1. The Group Sp(ri,]R). Let r(il) = fl + iJ he the tube over fl in Jc, which we identify with the 
space of complex n x n symmetric matrices. Let Go be the group of biholomorphic diffeomorphisms on 
T{n). Then Go is a Lie group with Lie algebra isomorphic to 5p(n,M) and acts homogeneously on T{fl). 
The group Sp(n, R) is isomorphic to a finite covering group of Go in the following precise way. The usual 
definition of Sp(n, M) is 

Sp(n,R)-{geSL(2n,R) IgjV 
where j = (^^-^ ^ SL(2n,M). Each element g € Sp(n,R) can be written in block form as 

A B 
C D 

where A, B, G, and D are real matrices. Defined in this way Sp(n,R) acts by linear fractional transfor- 
mation on the upper half plane J + iil. Let G be the group defined by 

This means then that the (1,2) and (2,1) entries of an element of G are purely imaginary matrices. 
Clearly G is isomorphic to Sp(n,R). It acts on the right half plane T[n) = + iJ hy linear fractional 

transformations: ii g — G G and z e T[fl) then 

g-z^ {Az + B){Cz + D)-\ 

It is also a finite covering group of Go. For g ~ e G we have the following relations among 

A, B, C, and D: 

A^C - G*A = AB^ - BA* = 

A*D - G*S = / AD* - BC* = I 

B*D - D^B = CD* ~ DC* = 

B*C-D*A = -I AD*~BC*^-I 

2.2. Some Subgroups of G. Let e be the n x n-identity matrix. Then e G i7 C T{fl). Let K be the 
stability subgroup of e in G. Then 

^^ = {(b ^) eGM±i?et/(n)|^f/(n). 

Here the isomorphism is given by 

'A B 



^B + 

The subgroup ii' is a maximal compact subgroup of G and G/K is naturally isomorphic to T{fl) by the 
map gK — > g • e. The connected component of the identity of the subgroup of G that leaves fl invariant 
is isomorphic to H via the map 



(h*)-' 



h 
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This map realizes L as a subgroup of G as well. In fact, we have 
via the above isomorphism. 

2.3. Lie Algebras. If P = then G = P-iSp(n,M)P. From this it follows that the Lie algebra, 

g, of G is given by g = P^^sp(n, K)P, and hence 

a lb 
—ic —a* 

We define a Cartan involution on g by 

e{x) ^ -X* . 

It induces a decomposition of g and gc into ±l-eigenspaces, g = 6+p and gc — tc+Pc- The +l-eigenspace 
t is the Lie algebra of K. These spaces are given by 



(2.1) 



e s[(2n,C) I a, 6,c real, 5 = 6*,c = c*| 



e = {X e g I 0{X)=X} = ||^"^ Gs[(2n,C) | a, 6 real, a = -a*, 6 = 6*, tr(a) = 
p = {X e g I 6'(X) = -X} = e 5l(2n,C) | a, 6 real, a = a*, 5 = 6*| . 



and 



Their complexifications are given by 

a b 
b a 



6c = -i ( ? ^ ) e sl{2n, C) \ a,b complex, a = -a^, b = 6^ tr(a) = 



and 



pc = <j ( ) "= s[(2?i, C) I a, 5 complex, a = ,b ^ b^ 



It is clear that iiTc acts on pc. This representation decomposes into two parts. For that let z = 
^ . Then z e 2(6c) and ad(z) : gc ^ gc has eigenvalues 0,2,-2. The 0-eigenspace is tc, the 
2-eigenspace is denoted by p^ and is given by 



a —a 
a —a 



e gc I a = a* )■ C Pc, 



and the —2-eigenspace is denoted by p+ and is given by 



p+ = <{ ( e gc I a = a* K PC 

Each of the spaces p^ are invariant under Kc and irreducible as representation. Note that this 
is not necessarily the standard notation. In our notation the eigenvectors (in p+) with — 2-eigenvalue 
correspond to annihilation operators while the eigenvectors (in p^) with + 2-eigenvalue correspond to 
creation operators. These operators will be described is Section 01 below. 

3. The Highest Weight Representations (tt^, 7Yy(r(rJ))) 

In this section we introduce the highest weight representations (tt^ , 7i^(r(r2))) and state the main results 
needed. We also introduce the Laplace transform as a special case of the restriction principle introduced 
in [El. 
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3.1. Unitary Representations of G in 0{T{Vl)). In this subsection we define a series of unitary 
representations of G on a Hilbert space of holomorphic functions on T{Q). These representations are well 
known. Let G be the universal covering group of G. Then G acts on T(0) by {g, z) i— » n{g) ■ z where 
K : G ^ G is the canonical projection. For v > n\ei T-iu{T{^)) be the space of holomorphic functions 
F : T{n) C such that 

(3.1) \\F\\l:^a^ [ \F{x + iy)\^A{xy-''''+^Uxdy<oo 

JTin) 

where 



2 

Then Ti.iy{T{Q)) is a non-trivial Hilbert space with inner product 

(3.2) {F\G)=a^ [ F{x + iy)G{x + iy) A{xY-^"+'^^ dxdy . 

For V < n this space reduces to {0}. If = n + 1 this is the Bergman space. The space Ti.i,(T{fl)) is a 
reproducing kernel Hilbert space. This means that point evaluation 

: n^Tin)) ^ c 

given by E^F = F{z) is continuous, for every z S TiVl). This implies the existence of a kernel function 
JC^ e n^{T{n)), such that F{z) = {F\K^) for all F g n^{T{n)) and z G T{n). Set K{z,w) = K^{z). 
Then K{z,w) is holomorphic in the first variable and antiholomorphic in the second variable. The 
function K{z,'w) is called the reproducing kernel for T-l^(T{n)). We note that the Hilbert space is 
completely determined by the function K{z,w). In particular, we have: 

(1) The space of finite hnear combinations Tiu{T{Q))° :— {J2'^j-^wj \ cj G C, wj e T{Q)} is dense 

in nATin)); 

(2) The inner product in Ti.,y{T{rt))° is given by 

^3^''"^ I '^^^^>' = Y CjdkK(zk,Wj) . 

j k J j,k 

We refer to [2 E] for more details. 

3.2. The Unitary Representations {7r^,n^{T{n))). For .g e G and z G T{n), let J{g,z) be the 
complex Jacobian determinant of the action of 5 G G on T(f2) at the point z. We will use the same 
notation for elements g G G. A straightforward calculation gives 

J(g,z) = det(Gz + D)-"-\ .9 = (^^ G G and z G T(f7). 

We also have the cocycle relation: 

J(ab, z) = J(a, b ■ z)J{b, z) 
for all a,b E G and z G T(f7). It is well known that for > n the formula 

(3.3) Mg)fiz) = J{g~\z)^^f{g~' ■ z) = dei{A - z Cyfig^' ■ z) 

defines a unitary irreducible representation of G. In [S], ^1], and it was shown that this unitary 
representation (tTj^, 7i^(r(ri))) has an analytic continuation to the half-interval v > {n — l)i. Here the 
representation ■Wy is given by the same formula H3.3|l but the formula for the norm in H3.1|l is no longer 
valid. There are also finitely many equidistant values of v that give rise to unitary representations, but 
they will not be of concern to us here. 
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In the following theorem we summarize what we have discussed and collect additional information 
from j2] and (cf. p. 260, in particular, Theorem XIII. 1.1 and Proposition XIII. 1.2). 

Theorem 3.1. Let the notation he as above. Assume that for v > n then the following hold: 

(1) The space Ti.i,(T{Q)) is a reproducing Hilbert space. 

(2) The reproducing kernel of TLv{T{^)) is given by 

(3) If V > i(7i— 1) then there exists a Hilbert space 'Hij{T{n)) of holomorphic functions on T (fl) such 
that Ki,{z,w) defined in (2) is the reproducing kernel of that Hilbert space. The representation 
TT^ defines a unitary representation of G on Ti.^. 

(4) If V > ^(ri — 1) then the functions 

A(z + e)->™ (^^-p^^ , me A, 

form an orthogonal basis of T-C,y{T{il))^ , the space of L -invariant functions in T-l^{T{Vl)). 

3.3. The Restriction Principle and the Laplace Transform. The restriction principle ^| is 
a general recipe to construct unitary maps between a reproducing kernel Hilbert space of holomorphic 
functions and L^-spaces on a totally real submanifold. Suppose M is a complex manifold and M{M) is a 
reproducing kernel Hilbert space of holomorphic functions on M with kernel K. Suppose X is a totally 
real submanifold of M and a measure space for some measure /x. 

Assume we have a holomorphic function D on M , such that D is positive on X, and such that the 

map 

R : H(M) -> L^{X, ^i) , 

given by Rf{x) = D{x)f{x), is densely defined. As each / is holomorphic, its restriction to X is injective. 
It follows that R is an injective map. We call R a restriction map. Assume R is closed and has dense 
range. If K{z,w) = Kw{z) is the reproducing kernel for EI(M), and / G L'^{X,d^), then 

R*f{z) = {R*f\K,) 

= U\RK.) 

= I f{x)D{x)K{z,x)dfi{x) . 
Jx 

In particular, if we set ^(x,?/) — D{y)D{x)K{y,x), then RR* is given by 

RR*f{y)^ f f{x)^{x,y)dii{x) 
Jx 

and thus is an integral operator. Consider the polar decomposition of the operator R* . We can write 

R* = UVRR*, 

where U is a unitary operator 

U : L^{X,^i) -> H(A/). 

The unitary map U is sometimes called the generalized Segal-Bargmann transform. In many applications 
of the restriction principle, M and X will be homogeneous spaces with a group H acting on both. When 
the restriction map R is _ff-intertwining so will the unitary operator U . This is exactly what happens in 
the situation at hand. Here we can take D = \ and define R : TLv{T{VL)) L^i^) by 

Rf{x) = fix) . 

Then we obtain the following: 
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Theorem 3.2. The map R is injective, densely defined and has dense range. The unitary part, U , of 
the polar decomposition of R* : R* — U\J RR* , is the Laplace transform given by 

Ufiz)=C,f{z)= I e-(^l^)/WdA'.(^) 

Furthermore, 



n 



C,{K{h)f)^TT,{h)C,{f) 

for all h E H . In particular, induces an isomorphism C^, : L^{il)^ 7Y^(T(il))^. Moreover 

C(C)=ro(m + i.)C. 

Proof. The first proof of this theorem was done for SL(2, M) in J^. The general case is on pages 187-190 
of El- □ 

Remark 3.3. Rossi and Vergne obtained the unitarity of the Laplace transform using a result of 
Nussbaum. 

The unitarity of the Laplace transform allows us to transfer the representation, tt,^, of G on Hu{T{ft)) 
to an equivalent representation of G on which extends Xi, by the above theorem, . We will denote 

the extension by Xi, as well. It is possible to describe on various subgroups of G whose product is 
dense in G. However, it is a difficult problem at best to describe a global realization of X„ on all of 
G. However, part of the point of this paper is to give a formula for the derived representation of Xi, on 
the Lie algebra of G and its complexification. It is from the derived representation that new differential 
recursion relations arise that relate the generalized Laguerre functions. 

4. The Action of qc 

In this section we introduce some subalgebras of sp(ri, C), the complexification of the Lie algebra of G, 
and explain how they act in the Hilbert space TLu{T{^)). 

4.1. The Derived Representation on 7ii,(r(il)). Denote by 'Hu{T{Vl))°" the space of functions F G 
7ii/(T(ri)) such that the map 

R 9 i i-^ 7rj.(exptX)i^ G n^{T{Vl)) 

is smooth for all X E Q ^ sp(n,M). If / e Cf{G), then 7r^(/)i^ = fig)7r^{g)F dg is in n^{T{n))°° 
and it follows, that 'H„{T{il.))°° is dense in Ti^{T{^)). The Lie algebra representation, denoted also by 
TTu, of g on 'Hi,{T{n))°° is given, by differentiation as follows: 

t^o t 

= ^7r,(exp(tX))F|t=o- 
at 

Note that the limit is taken in the Hilbert space norm in 7i^(T(ri)), but it is easy to see that if F e 
liy{T{n))°°, then in fact for X E q: 

(4.1) ^AX)F{z) = ^ J(exp(-tX), z)^iFieM-tX) ■ z)\t=o , 

for all z e T{n). We extend this by complex linearity to gc- 
Define D^j by 

D^F{z) = ^F{z + tw)\t^Q = F'{z)w 
dt 

where F' denotes the derivative of F. 
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Lemma 4.1. Suppose z,w are n x n matrices over C and z is invertible. Then 

D^detiz)"- = 7idet(z)"tr(z~^w). 
Proof. This follows from the chain rule and the fact that 

det(z) = det(z + tw)\t=o = det z-^ det(l + tz~^w)\t=o = det(z) tr^z^^w) . 



□ 



The following proposition expresses the relevant formulas on ic, p+, and p . It's proof is a straight- 
forward calculation using Lemma l4.ll 



Proposition 4.2. For each piece of the Lie algebra of qc introduced in Suhsection \2.t\ we have 
(1) MX)F{z) = vtr(bz)F{z) + D,a-az-b+zb.F{z), X^(l '!) eic 



(2) 7r,{X)F{z) = -vtr(az + a)F{z) - D(,a+az)+zaz+aF{z), X = (^^^ € p+ . 

(3) Tr,{X)F{z) = -i/tr(-az + a))F{z) + D_(^,,+az)+zaz+aF{z), X ^ {^^ G 

Remark 4.3. We note that these equations are closely related to the Jordan algebra structure of J . In 
particular, we have that za + az = 2z o a, where o denotes the Jordan algebra product aob = ^{ab + ba). 
Furthermore zaz — P{z)a where P denotes the quadric representation (cf. p 32). Finally we note, 
that tv[az) — tr(z o a). These formulas therefore clearly indicate that the more general results are 
expressible in terms of Jordan algebraic constructs. 

4.2. Highest Weight Representations. The fact that tt^ is a highest weight representation plays a 
decisive role in the recursion relations that we obtain. At this point we explain what this notion means. 

We assume G is a Hermitian group, which means that G is simple and the maximal compact subgroup 
K has a one dimensional center. The Hermitian groups have been classified in terms of their Lie algebras. 
They are su(p, q), sp(n,M), so*(2n), 5o(2,n), and two exceptional Lie algebras. The assumption that K 
has a one dimensional center implies that G/K is a bounded symmetric domain. In particular G/K is 
complex. It also implies that the complexification of the Lie algebra, gcj has a decomposition of the form 
gc = P"*" ffi ffi P~, Specifically, p+, 6c, and p~ are the —2, 0, 2-eigenspaces of ad(z), respectively, where 
z is in the center of 6c- 

Lemma 4.4. We have the following inclusions: 

[«c,P*] cp± 



Proof. Suppose X ep^,Y £p~ , and Z etc- Then 

ad(z)[X,r] = [ad(z)X,y] + [X,ad{z)Y] = -2[X, + 2[X, F] = 0. 

This implies [X, Y]etc- Similarly, 

ad(z)[Z,X] = [ad(z)Z,X] + [Z,ad(z)X] = -2[Z,X]. 

This implies [Z, X] e p^. A similar argument shows that [Z, F] e p^. □ 

Suppose TT that is an irreducible representation of G on a Hilbert Space H. We say tt is a highest 
weight representation if there is a nonzero vector w S EI such that 

t:{X)v = 0, 
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for all X e p+. Let Ho be the set of all such vectors. The following theorem is well known. 

Theorem 4.5. Suppose tt is an irreducible unitary highest weight representation of G on M and Ho is 
defined as above. Then (7r|;<-, Ho) is irreducible. Furthermore, there is a scalar X such that 

7r(z)i> = Xv, 

for alive Ho. // 

H„ = {w e H I n{z)v = (A + 2n)v} , 

then 

H = 0H„. 

n>0 

Additionally, 

^(Z):H„^H„, Zetc; 
^(X) :H„^H„_i, Xep+; 

^(y) :H„ ^H„+i, yep-, 

where, in the case n = 0, H_i is understood to be the {0} space. 

Proof. By Lemma |4.4I Ho is an invariant K-space. Suppose Vo is a nonzero invariant subspace of Ho 
and Wo is its orthogonal complement in Ho. Define V„ inductively as follows: 

V„ = span{^(y)w \Ye\>-,ve Yn-i} ■ 

Let V = eV„. Define W„ in the same way as V„ and let W = ®W„. Then, by Lemma V and 
W are invariant qc subspaces of H. Since tt is unitary V and W are orthogonal. However, since tt is 
irreducible and V is nonzero, it follows that V = H and hence W = 0. This implies Wo = and thus 
ttIk is irreducible. Since 7r(z) commutes with tt{K) Schur's lemma implies that 7r(z) = A on Ho for some 
scalar A. Since Vo = Ho, induction. Lemma 14.41 and irreducibility of tt implies that V„ = H„. The 
remaining claims follow from Lemma 14.41 □ 

Remark 4.6. The operators 7r(X), X e p"*", are called annihilation operators because, for v in the 
algebraic direct sum 0H„, sufficiently many applications of tt{X) annihilates v. For Y E the 
operators 7r(y) are called creation operators. 

Remark 4.7. A straightforward calculation gives 

for all X £ and that Hy(T(f2))o = Cqq. Thus (7r,y, H^(T(f]))) is an irreducible unitary highest weight 
representation of G and by unitary equivalence so is (A^, L^(f2)). 

5. The Realization of X^ Acting on L'^{il,d^j,^) 

In this section we determine explicitly the action of gc on L^{^1, d^v). More specifically, we define Xi, via 
the Laplace transform by the following formula 

x,{x)^cz^^,Ax)c, 

and will determine exphcit formulas for A^(A), for X e p+, X e 6c, and X e p^. 



LAGUERRE FUNCTIONS ON SYMMETRIC CONES 



13 



5.1. Preliminaries. Let Eij be the n x n matrix with a 1 in the position and O's elsewhere. 

Define Eij = ^{Eij + Ej^i). Then the collection ^Eij | 1 < i < j < n| is a basis of J and Jc, the 

real and complex symmetric matrices. Furthermore, (^Eij\Ek,i^ — ^{SjkSu + SjiSik), which implies 

^Ei j I 1 < j < J < is an orthogonal basis. Set Dij = Dj^, . and observe that Dij — Dj^i. The 
gradient of /, V/, is defined by 

{y.f{x)\u)^DJ{x). 

Proposition 5.1. Suppose f,g^ L^{Q.,diJLu) are smooth and f vanishes on the boundary of the cone fl. 
Let 1 < i, j < n. Then 



(1) 



Dt,jf{s)g{s)ds ^- f{s)Dijg{s)ds 



(2) 

' e-(^l^)z,.,/(s)ds = / e-(^l^)A.,/(s)rfs. 
! Jn 

Proof. (1) is Stokes Theorem and (2) follows from (1) and the fact that Di^e^^^''*^ = — e~^^l*-'zi j, 

z e Jc. □ 

5.2. The Representation A^. Recall that we determined the action of tp, p+ and on 7Y^(T(f7))°° 
in Proposition We denote the subspace of smooth vectors in Ll{Q) by Thus / e Ll{ri)°° 

if and only if the map 

R3t>-^ K{exptX)f e Ll{n) 

is smooth for all X G g. Thus 

The action of g on Ll{n)°° is, as usual, defined by 

A,(X)/ = lnn , 

for X e g, and then by complex linearity the action extends to gc- The following theorem collects the 
corresponding equivalent action on the Hilbert space We remark again that these formulas can 

be stated in terms of the Jordan algebra structure of J indicating the extension of these results to other 
tube domains. 

Theorem 5.2. For f e L't{n) a smooth function we have: 

(1) X^{X)f{x) = tr[{bx + (ax - xa ~ vb)\J ~ x\Jb\J]f{x), X ^ a) ^ '^^^ 

(2) K{X)f{x) = tv[{va + ax+ {ax + xa + ua)W + xWaW]f{x), X = e p+ 

(3) X,,{X)f{x) = tr[(j/a - ax + [ax + xa - i^a)\/ - x\/a\/]f{x), X ^ T Z^) ^ 
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5.3. The Case of 5p(l,R). The proof, which appears in the appendix, is very long and computational. 
However, to convey the main ideas of the proof we will discuss the simpler case of Sp(l,M), which is 
isomorphic to SL(2,R). A detailed account of this case (modelled on the upper half plane) is found in 

m 

Let G = I I e SL(2,M)| . The group G acts on the right half plane T(R+) by 

linear fractional transformations. The complexification, gc, of the Lie algebra of G is sZ(2,C), a three 
dimensional Lie algebra spanned by 



and Y 



1 -1 
1 -1 



Proposition 14 . 21 for this case reads: 

Theorem 5.3. The action o/s[(2,M) on the right half-plane is given by: 

(1) 7r^(z)i^(z) = i^zF{z) + (z2 _ l)F'(z) 

(2) n,{x)F{z) = -i^{z + l)F{z) - {z + l fF'{z) 

(3) My)F{z) = iy{z - l)F(z) + (z - ifF'iz) 

To find the corresponding action on L^(M+) we must compute the operators that corresponds to D^, 
M^, M^2, o and M^i o in L^(M+)°°. Here M stands for "muhiplication operator". To do this, 
requires several uses of integration by parts, a special case of Stokes theorem. It was exactly this kind of 
computation that was done in ^ and we repeat it here: 

For Dz we have: 



Thus < — > M_f 
For Mz we have 



dz Jq dz 

= U-tf{t)) 



zUf){z) = / -^-f{t)t-^dt 
e-''^^ifitr-')dt 



Thus Mz < — > D + Mu^ . 

t 

We calculate similarly and get 



< y + ^— '-D 

Thus: 



) 1 

e-'-\f'it)+'^f{t))f^-'dt. 



t^ 



Lemma 5.4. Let the notation be as above. Then the following holds: 
(1) DzoC,^-C,o{Mt); 



(2) MzoC^^C^o{D 

(3) Mz oDzoL^ = C^o {-tD - v); 

(4) Mz2 oDzO L^^ L^o (-tD^ - 2vD 



t 



Combining Theorem 15 . 31 and Lemma l5 . 41 gives 
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Lemma 5.5. Let the notation he as above. Then the following holds: 

(1) A^(z) = -tD^ -lyD + t; 

(2) A^(x) = tD^ + + 2t)D + + t); 

(3) A^(y) = -tD"^ + {-V + 2t)D + {iy-t). 

Note that this is Theorem 15. 21 for this special case. 

One more ingredient is necessary for determining the classical recursion relations. This is a direct 
calculation and given in the following lemma. We note at this point that such a direct calculation is not 
done in the general case; deeper properties of the representation theory must be used. (c.f. Proposition 

O) 

Lemma 5.6. Let q'^^{z) ^ {z + l)-" [j^] ■ Then 

(1) 7r,(z)g^, = (z. + 2m)C 

(2) 7r^(x)g^ = -2mq'^__^ 

(3) 7r,(Y)C = 2(z^ + m)C+i 

The combination of Lemmas 15.51 and 15.61 gives the classical recursion relations stated in the intro- 
duction and proves Theorem IS. 31 for the classical case. 

6. Differential Recursion Relations for 

We now turn our attention to differential recursion relations that exist among the generalized Laguerre 
functions. These relations are obtained by way of the highest weight representation and generalize 
the classical case mentioned in the introduction. 

We begin with some preliminaries and a result found in [2]. First we notice that in general the Lie 
algebra gc does not map {L'l{U)°°)^ into itself. For the Laguerre functions the full Lie algebra is too big; 
we will in fact only need the much smaller Lie algebra g^, which maps {Ll{^)°°)^ into itself. It is well 
known, that in case g is simple, then Q^i ~ s[(2, C). We choose z, x, Y so that the isomorphism, which we 
will denote by (p, is given by 

! J) ' (-1 -i)' (1 ~\ 

Furthermore, we can assume that </3(X*) ~ ip{XY . This shows that several calculations can in fact be 
reduced directly to s[(2, C). We will come back to that later. 

Define := i(x + y). Then ip{Z°) = and Z° is in the center of f). For m G A let 

^ (^-^ ^ TT ^3 -^k- \{i + 1 - fc) 
iX - mk - \{j - k) 

Then by Lemma 5.5 in 2 we have: 

Proposition 6.1. The action of z and Z° is given by: 

(1) 7r,(z)g^ = (ni. + 2|m|)C- 

(2) M-2Z")q-^ = (m™^>m-e, - Y^Ui^ + - \{j - l)))Cm(j)C+e, • 

Corollary 6.2. Let the notation he as above. Then the following holds: 

(1) A.(z)C = (ri;^ + 2|m|)C. 

(2) A.(-2Z°)C = E;=1 (™'"e,)K + l))C-e, - E;=1 C™(j)C+e, • 
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Proof. This statement follows from Proposition 16.11 and the following facts: Xu{X) = £'^^tt,^{X)C^, 
Cy{l'^ — ro(m + v)q'^, and Proposition ^21 In each of these formulas if either index m + or m — 
is not in A then it should be understood that the corresponding function does not appear. □ 

Theorem 6.3. The Laguerre functions are related by the following differential recursion relations: 

(1) tr(-a;VV - i^V + xyi^{x) = {nv + 2|m|)C(a;)- 

(2) tr(xVV + {vl + 2x)V + {vl + x))l-^{x) = -2E;=i (^"'ej) K - 1 + - (j " l))C-e,(^) 

(3) tr(-a;VV + {-vl + 2x)V + {vl - x))r^{x) = 2 E^^i Cm(j)C+ej {^)- 

Proof. If z = then substituting a = and b = \ into Proposition 15. 21 part 1, gives 

A^(z) = tr(-a;VV - + x). 

Combining this with part 1 of CoroUarv 16 . 21 gives the first formula. 

Recall that x e p+ and Y e p^. According to Theorem 14. 51 we have that Ajy(x)£J^ has to be a linear 
combination of , with m'^ < ruj for all j. Similarly, has to be linear combination of those P^, 

with m'j > rrij. The statement follows now from CoroUarv 16.21 and the fact that 2Z° = x + Y. □ 



7. Some Open Problems 

There are still several open question that require further work. We mention three of these. One is 
a relation to the classical Laguerre polynomials, the other two are natural generalization of classical 
relations. 



7.1. Relation to Classical Laguerre Functions. Every positive symmetric matrix A can be written 
ds A — kDk^^, where k G SO(n) and D — d{ti, . . . , t„) is a diagonal matrix with tj > 0. Thus, if 

ill = {d{t) I t e (R+)"} ~ (R+)" 

then 

Lfli. 

As the Laguerre functions are i-invariant, it follows that they are uniquely determined by their restriction 
to Qi. Let T{ni) := {d(x) + id(y) | x e (M+)",y e M"}. Then r(rJi) ~ (R+ + and the group 

SL(2,IR)" acts transitively on the right hand side. But it is well known, that SL(2,IR)" can be realized 
as a closed subgroup of Sp(n,M). It follows therefore, that the generalized Laguerre functions can be 
written as a finite linear combinations of products of classical Laguerre functions. It is a natural problem 
to derive an exact formula. 



7.2. Relations in the A-parameter. It is well known that the classical Laguerre polynomials satisfy 
the following relations: 

xL^ = (n + A + l)i^i-(n+l)i^-i 
xL^^ - (n + A)i^li-(n-a:)i^i 



In ^ it was shown, that these relations follows directly from the representation theory of s[(2,M). It is 
therefore natural to look for similar relations for the generalized Laguerre polynomials and functions. 
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7.3. Relations in the x, y Parameters. Several other classical relations should be extended to the 
general case. We name here only the following 

m 

This relation is closely related to the decomposition of the tensor product of two highest weight represen- 
tations and we expect that a similar relation can be derived also for the general case. Notice, however, 
that for general Laguerre polynomials the right hand side is L-invariant in the x and y variable while 
that is not the case on the left hand side. Thus any generalization will involve a projection (averaging 
over L) onto the L-invariant functions. 
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8. Appendix: Proof of Theorem (|5.2|1 



Proof. We will prove the theorem for the case x G p . The other two cases are done similarly. For 

rt+l 
2 ■ 



convenience we let to = 



a —a 



Let x = ( ) e p+. By Proposition l4.2l 



Let 

A = -Daz+zaCyf{z) 

B = -DaC,f{z) 

C = -D,azC,.f(z) 

Calculation of A. 

A = -Daz+zaCuf{z) 

= - [ Daz+zae'^'-^^^f{x)det{xr dx 

= [ e-'^'-^'='>{az + za\x)f{x)det{x)'^ dx 
Jn 

= / e~''''^''^aikZkj + Zikakj){xjif{x)det{x)"^) dx 

= y. I e-("l")(a.feDfe, + au,D,u){x,,f {x)det{xD dx 

= V / e~(^l=")(ajfc-(4i'5ji + htSjj) + akj-{SkjS^^ + Sij6ki))f{x) det(a;)™ dx 
+ e^''^^''H(^^kXj^Dkj f (x) + akjXj^D^k fix)) det{x)"^ dx 

. / e'^^'''^'-'^f(x)mdei{x)"'{aikXj^ii{x^^Ekj) + akjXj^ii{x'''^Eik))dx 



(n + l)tr(a) / g-^^'^VCa;) det(x)'" 



e-(^l-)((ax + xa),kD,kf{x)) det(x)'" dx 

! 

+ / e-(^I^V(a;)TOdet(x)'"tr(x-i(xa + ax))dx 

= (2to + n + 1) tr(a) [ e''^''^'=^ f{x) det{x)"' dx + [ e-(^l^)(tr((ax + xa)V/)(x)) det(x)™ dx 

= 2i/tr(a) / e-(^l=^V(a;)det(x)" dx+ / e~(^l^)(tr((ax + xa)V/)(x)) det(x)'" dx 
Jn Jn 

= 2vir{a)L^{f){x) + C^{{tr{ax + xa)V)/)(x) 
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Calculation of B:. 

B = -DaCufiz) 

= [ e-(^l^)(ak)/(a;)ciet(a;)™da; 
Jn 

= [ e-(^I^Ur(aa;)/(a;) det(a;)™rfa; 
Jn 

= jCt,{tT{ax)f){x) 

Calculation of C:. 

C = -D,azC.f{z) 

= I e-^''\''\zaz\x)f{x)Aei{x)"'dx 
Jn 

= ^ / e~^''^''^akiZikZijXjif{x)det{x)'^dx 

= E [ e-^'^''^akiZikDij{xjif{x)det{xr)dx 

= E / e-^'^^^akiZikl-{5ij5ji + Sii5jj)f{x)det{xrdx 

+ I e-('^^^ak^z,kXJ^{D^Jf){x)det{x^dx 

+ E / akiZikXjif {x)m det{x)"^ tv{x~^ Eij)dx 

= E / e-^'\='^akizik{'^6u)f{x)detixrdx 

+ E / e-^'\''^akiDik{xjiDijf{x)det{xDdx 

+ E / e~'^''^''^akiZikXjif{x)mdet{x)'^ti{x~^Eij)dx 

= ^^tr(az) / e-(^I^V(a;)det(x)™dx 

+ E / e-^'\^^akihSij5ki+Sii5kj)iDijf){x)detixrdx 

+ E / e-^'^^^akiXji{DikDijf){x)det{xrdx 

+ E [ e-^'\='^akiXjiDijf{x)mdet{xrtv{x-^Eik)dx 
i,3,ky^ 

+mV / e~'^''^^'^aikZkiSiif{x)det{x)'^dx 
j,k 
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2 in 

+ V / e-^'\='^akA{Sij6ki + 6u6kj){Dijf){x)det{xrdx 

+ V / e-(^l^)awa;,i(AfcA,/)(a;)det(x)'"da; 

+ f e-^^^^^akiXjiDijf{x)mdet{xrtT{x-'Eik)dx 

+TnV / e~^''^''^aikZkiduf{x)det{x)"'dx 

= vtviaz) f e-(^I^V(a;)det(a;)'"rfa; 

/ e-(^l^)(tr(aV)/)(a;)det(a;)'"da; 

+ / e-(^l^)((tr(a;VaV)/)(a;)det(a;)'"da; 
Jn 

= i'ti{az)C^f{x) + ujC„{ti{aV)f){x) + jC^{tr{xVaV)f){x) 



■Ku{x)Cvf{z) = -utT{aZ + a)C^f{z) - Daz+a+zaz+zat^vf{z), X € p + 

= -vtT{az + a)C„f{z) + 2iytT{a)C^{f){x) 

+C,{{tr{ax + xa)V)f){x) + £,(tr(ax)/)(x) 

+i^tT{az)C^f{x) + i^C^{tv{ay)f){x) 

+£^(tr(xVaV)/)(a;) 
= iytT{a)C4f){x) + jC4{tr{ax + xa)V)f){x) 

+C4tv{ax)f){x) + ujC^{tr{aV)f){x) 

+C,{tr{xVaV)f){x) 
= £,/(tr(;ya + ax + {ax + xa + va)\/ + xVaV)/(a;)) 
Taking the inverse Laplace transform of each side gives the desired result. 

□ 
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